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1 Introduction 

The study of mesoscopic phenomena is one of the most active areas of today's 
solid state physics. One can observe signatures of mesoscopics in a large num- 
ber of different physical systems, and a comprehensive review would not be 
appropriate in the present context. The present volume focuses on Kadanoff- 
Baym-Keldysh Green functions, and consequently wc restrict ourselves to cer- 
tain examples which have been studied with the help of these techniques. The 
generic system we have in mind is a semiconductor heterostructure where 
charge carriers are introduced either by modulation doping, or they flow in 
and out of the system through metallic or superconducting contacts. Transport 
physics in these systems can roughly be divided into two categories: perpendic- 
ular transport and parallel transport, according to whether the charge carriers' 
motion is perpendicular or parallel to the layers that form the heterostructures. 
A representative example of perpendicular transport is the resonant-tunneling 
diode (RTD), which consists of alternate layers of semiconductor materials with 
different band gaps; a schematic conduction band diagram is shown in Fig.|l]. 
Charge carriers entering from left may, at a certain bias voltage, be tuned to 
the quasibound state in the quantum well, which results in a large enhance- 
ment of the transmitted current. At off-resonance conditions only a small 
current can flow, because transmission through the classically inaccessible re- 
gions is exponentially suppressed. This leads to a nonmonotic current-voltage 
characteristic, and a number of device applications have been proposed, whose 
operating principles are based on this property. 

In the case of parallel transport much attention has been devoted to quan- 
tum point contacts (QPC), see Fig.||for a typical experimental configuration, 
and other structures based on similar ideas. Here the key ingredient is metal- 
lic gates that are deposited on the heterostructure; by adjusting the gate po- 
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Figure 1: Double- barrier semiconductor heterostructure biased close to resonance, where 
charge carriers emerging from the source contact are matched to the energy of the quasibound 
state £0 in the quantum well. Occupied contact states are shown as hatched, and the band 
bending is due to charge accumulation or depletion. 




Figure 2: Coulomb island, which consists of two tunable quantum point contacts QPCl 
and QPC2, and a side gate which allows one to vary the chemical potential, and hence the 
charge density in the central region. The two-dimensional electron gas underlying the gate 
structure is depleted outside the hatched regions. 

tentials it is possible to deplete the underlying two-dimensional electron gas, 
and thus introduce spatial modulations of the two-dimensional charge density. 
Quantum point contacts are based on a split gate geometry: here, at sufh- 
ciently high negative gate voltages, the effective connection between the two 
unmodulated electron gases ( "source" and "drain" ) is so narrow that perpen- 
dicular mode quantization becomes observable, and the measured conductance 
is an integer multiple of the quantum unit of conductance, /h. In later sec- 
tions we describe simple models pertaining to structures like the one shown in 
Fig.|. 

The hallmark of mesoscopic phenomena is the phase coherence of the 
charge carriers, which is maintained over a significant part of the transport 



process. The interference effects resulting from this phase coherence are re- 
flected in a number of experimentally measurable properties. Weak localization 
can be understood as an increased return probability (and hence increased re- 
sistance) due to coherent backscattering of charge carriers. Another example 
where phase coherence is central is the Aharonov-Bohm effect, where interfer- 
ence of two different transport paths in a ring geometry results in an oscillatory 
magnetoresistance. Yet another example is universal conductance fluctuations, 
where the conductance of a sample displays rapid changes on the scale of e^/ft. 
(hence the "universality" ) when an external control parameter is changed. The 
external parameter could be magnetic field, or thermal cycling, and the fluctu- 
ations reflect changes in the conducting channels either due to different impu- 
rity configurations (thermal cycling), or differences in the way the conduction 
channels are located in the sample (magnetic field). In these notes we focus on 
an alternative way of affecting the phase coherence: external time- dependent 
perturbations. The interplay of external time dependence and phase coherence 
can be phenomenologically understood as follows. If the single-particle ener- 
gies acquire a time dependence, then the wavefunctions have an extra phase 
factor, ip ~ exp[— i / dt's{t')]. For a uniform system such an overall phase 
factor is of no consequence. However, if the external time dependence is differ- 
ent in different parts of the system, and the particles can move between these 
regions (without being "dephased" by inelastic collisions), the phase difference 
becomes important. 

The interest in timc-dcpcndcnt mcsoscopic phenomena stems from recent 
progress in several experimental techniques. Time dependence is a central in- 
gredient in many different experiments, of which we mention the following: (i) 
Single- electron pumps and turnstiles. Here, time-modified gate signals move 
electrons one by one through a quantum dot, leading to a current which is 
proportional to the frequency of the external signal. These structures have 
considerable importance as current standards. The Coulombic repulsion of 
the carriers in the central region is crucial to the operational principle of these 
devices, and underlines the fact that extra care must be paid to interactions 
when considering timc-dcpcndcnt transport in mcsoscopic systems, (ii) ac 
response and transients in resonant-tunneling devices. Resonant tunneling de- 
vices (RTD) have a number of applications as high-frequency amplifiers or 
detectors. For the device engineer, a natural approach would be to model 
these circuit elements with resistors, capacitances, and inductors. The ques- 
tion then arises as to what, if any, are the appropriate "quantum" capacitances 
and inductances one should ascribe to these devices. Answering this question 
requires the use of timc-dcpcndcnt quantum-transport theory. 

A central issue will be the treatment of interactions in the mcsoscopic 



3 



region, and, as we shall see, nonequilibrium Green function techniques are well 
suited for this purpose. First analyses of tunneling problems in nonequilibxium 
systems were performed already in the 70s by Caroli and co-workers Boa'cl. In 
later years the steady state situation has been addressed by a large number 
of papers, however the literature on time-dependent nonequilibrium transport 
treated with Green functions is much more restricted than in the stationary 
case (see the recent text-book by Haug and JauhoQ for extensive references for 
works published before 1996; this review will concentrate on some of the more 
recent advances. Additional useful and related information can also be found 
in the monograph by Dattaj.). 

Among the central results obtained with the nonequilibrium Green func- 
tion techniques is that under certaia|-conditions (to be discussed below) a 
Landauer-type conductance formula ltQ can be derived. The Landauer for- 
mula relates the conductance g of a mesoscopic sample [which is connected 
via "ideal" leads to two (or more) reservoirs] to its transmission properties, 
g = [e^ /h)T, where T is the quantum mechanical transmission coefficient of 
the sample. Conductance formulae have played an important role in the anal- 
ysis of many mesoscopic transport phenomena, and it is therefore of interest 
to investigate whether interactions and/or time-dependence can be treated in 
a similar fashion. This study forms the core of the present review. 

2 Nonequilibrium formulation of tunneling physics 

The total current measured in an external ammeter can be split into two 
contributions: the current flowing into the mesoscopic region, and the current 
flowing in and out of the accumulation or depletion regions in the leads, i.e., the 
displacement or capacitive current. We begin by deriving an expression for the 
tunneling current, and comment on models for the displacement current later. 
One should also note that the displacement current does not contribute to the 
time averaged current, which is often the experimentally relevant quantity. We 
shall also assume that the electric fields in the leads are effectively screened out, 
so that the voltage drop occurs across the mesoscopic region. This assumption 
sets an upper bound for the external driving frequency, which should not exceed 
the plasma frequency. Estimates for the cut-off frequency are in the range of 10 
GHz - 1 THz, depending on the device geometry and dimensionality B, which 
are sufficiently high for most present applications. 

We recall that the basic difference between construction of equilibrium 
and nonequilibrium perturbation schemes is that in nonequilibrium one can- 
not assume that the system returns to its ground state (or a thermodynamic 
equilibrium state at finite temperatures) as t ^ -|-oo. Irreversible effects break 
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the symmetry between t = —oo and t = +00, and this symmetry is .^''C'.'^Y.'V 
exploited in the derivation of the equihbrium perturbation expansion E2IEiHi3. 
In nonequihbrium situations one can circumvent this problem by allowing the 
system to evolve from —00 to the moment of interest (for definiteness, let us 
call this instant tg), and then continues the time evolvement from t = to back 
to i = — 00. The advantage of this procedure is that all expectation values are 
defined with respect to a well-defined state, i.e., the state in which the system 
was prepared in the remote past. The price is that one must treat the two 
time branches on an equal footing. 

In the context of tunneling problems the nonequilibrium formalism works 
as follows. In the remote past the contacts (i.e., the left and right lead) and the 
central region are decoupled, and each region is in thermal equilibrium. The 
equilibrium distribution functions for the three regions are characterized by 
their respective chemical potentials; these do not have to coincide nor are the 
differences between the chemical potentials necessarily small. The couplings 
between the different regions are then established and treated as perturbations 
via the standard techniques of perturbation theory, albeit on the two-branch 
time contour. It is important to notice that the couplings do not have to be 
small, e.g., with respect level to spacings or k^T, and typically must be treated 
to all orders. 



3 Current formulas 

We can now present the mathematical formulatioixof the problem. A detailed 
presentation can be found in a recent text-book El, and here we present just 
some of the central ideas. The contacts are assumed to be noninteracting, and 
the single-particle energy in lead a is given by 

£k,a{t) = ek^a + Aa(i) , (1) 

where Aa{t) is the external time modulation. The leads are connected to 
the central (or, mesoscopic) region via a hopping term with matrix element 
Vka;n{t), whcrc n labels the eigenstates of the central region. Collecting the 
various terms results in the Hamiltonian H = Hl + Hr + Ht + Hcen, or, 
explicitly: 

H = ^ek,ait)cl^^Ck,a + ^ [VfcQ;„(i)c^d„ + h.C. + iJcon [{^n}, {dJi}, i] , 

k,a k,a;n 

(2) 

where the central part Hamiltonian must be chosen according to the system 
under consideration. The operators {dn}, {dj^} refer to a complete set of single- 
particle states of the central region. The derivation of the basic formula for 
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the time-dependent current does not require an explicit form for i?con; the 
actual evaluation of the formula of course requires this information. We write 
^^cen = X^ra ^n{t)dl^dn + Hint, whcrc iJint could bc clcctron-phonon interaction, 
or an Anderson impurity: 

na q 

i?4 = C/E<T^".T4arfn,i- (4) 

n 

According to the basic ideas of the tunneling approach presented above, the 
occupations of the leads are determined by equilibrium distribution functions. 
Thus the Green functions for the contacts are known explicitly: 

= */(£L) exp [ - i y dtiEkcih)] , (5) 

= Tie{±t T i') exp [ - z y^ dtiefca(ii)] ■ (6) 

We start the derivation by considering the current leaving the, e.g., left contact, 
and entering the central region: 

JlW = (/lW) = {{-e))NL(t)) = -ie{[H,NL]) . (7) 

The commutator [H, Nl] is readily evaluated, and one finds 

Mt) - ^Re{ E ^fc",„(t)G<,Ji,<)} , (8) 

k,a,n 

which involves the time-diagonal part of the correlation function 

GtkJt^t') ^ z{cljt')d,,{t)) . (9) 

The next step consist of writing down the equation-of-motion for the time- 
ordered function ^.^{t, t'), and a subsequent analytic continuation with the 
Langreth rules 111 leads to 

(10) 



Substituting the expressions for the contact Green functions, Eqs.(||)-(||) fi- 
naUy yields 

JlW - -| f dt,j ^lmTr{e-(*-*)r^(6,ti,i) 

X [G<(t,ii) + /0(e)G'-(t,ii)] } . (11) 

Here the Green functions G^^'' are matrices in the indices (m,n), and the 
hnewidth functions T are defined as 

rti 

[r^(e,ii,i)]^„ = 2^^p„(e)K,„(e,i)F„V(e,ii)exp[-z / dt2Aa(e,t2)] , 

. ^^^^ 

where ^^(e) is the density of states. Equation ( |ll| ) is the main formal result 
of this report: the subsequent sections are devoted to exploring its special ap- 
plications. It is important to note that the current formula only involves the 
Green function of the central region. However, G^{t, ti) must be calculated in 
the presence of the coupling to the leads, which is a highly nontrivial task for 
an interacting system. Thus, Eq.(|ll|) can be viewed as a rather formal state- 
ment, but nevertheless it provides under many circumstances a very convenient 
starting-point for further calculations. 

4 Special cases 

4-1 Stationary limit 

In the stationary liiniJ. Eq. ( p^ ) can be further simplified and we get the result 
first reported in Refo: 

+ [/^(e)r^(6) - fnie)T%)] [G^e) - G''(e)] } (13) 
= jI J[/L(£)-M£)]T(e), (14) 

where 

Equation (|l^) holds for the special case when the couplings between the left and 
right leads are proportional (a constant coupling, occurring in the wide-band 
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limit is a frequently encountered special case) . The ratio in the curly brackets 
in Eq.(15) is well-defined because the F- matrices are, by assumption, propor- 
tional. The difference between the retarded and advanced Green functions is 
essentially the density of states. Despite the apparent similarity of Eq.jT^) to 
the Landauer formula, it is important to bear in mind that, in general, there 
is no immediate connection between the quantity T(e) and the transmission 
coefficient r(e). In particular, when inelastic scattering is present, there is 
no such connection. These expressions for the stationary current have a wide 
range of applicability for calculations of current- voltage relations in mesoscopic 

euctures. Recent applications inclujdc nonlinear currei^val±ape|Calculations 
, transport in carbon nanotubesll3'EJ, Kondo physics ll3tJ' — 'Eii23, Andreev 
scattering in semiconductor-superconductor hybrid systems (we return to 
these two topics below), tunneling 
of excitons in disordered systems t 



quantum dots or artificial moleculesE 



rrough magnetic barriersEa, delocalization 
aQc^gj^^^ions in transport thr^ 



1, physics of.nanpwiresi 



analysis of STM (Scanning Tunneling Microscope) experimentsESta, and many 
others. 



4-2 Noninteracting electrons 

Now iJcen = eo{t)(i^d = [eo + /S.Q{t)]Sd (we examine only the single-level case; 
a generalization to many levels and/or spin is straightforward, all the results 
are formally the same but the Green functions, self-energies etc. must be 
interpreted as matrices), and we can give the following jexplicit results for the 
occupation N{i) = —iG^{t,t), and the current Ji^/ji{t)U: 



L/R 

JL/R{t) = -^rL/R\N{t) + - I de/0/^(6)Im{Ai/^(6,t)} 



(16) 
(17) 



where 



^L/Rie.t) 



k— — oo 



eo-kj + iT/2 



(18) 



where Jk is the k:th order Bessel function, and we assumed a harmonic time- 
variation: Aq (t) — cos (7t) . Figures I and I give a numerical example. We 
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Figure 3; as a function of time for harmonic modulation for a symmetric structure, 

= Fji = r/2. The unit for the time-axis is h/F, and all energies are measured in units of 
r, with the values /il = 10, /ir = 0, eo = 5, A = 5, Al = 10, and = 5. The modulation 
frequency is uj = 2T /H. 



0.1 




0.0 



12 3 4 5 6 
Time 

Figure 4: The time-dependent current J{t) for harmonic modulation. The dc bias is defined 
via /ii = 10 and fin = 0, respectively. The dotted line shows (not drawn to scale) the time 
dependence of the drive signal. The temperature is kgT = O.IF. 



draw attention to the maxima in the plot for | Ap; these are related to photonic 
side-bands occurring at e = eg i few ED. More work afong similar lines can be 
found, e.g., in ReftS Bearing in mind the complex structure seen in Fig. |^ 
it is not surprising that the current in Fig. ^ displays a non-adiabatic time- 
dependence. The basic physical mechanism underlying the secondary maxima 
and minima in the current is the line-up of a photon-assisted resonant tunneling 
peak with the contact chemical potentials. 
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4-3 Average current 

For constant couplings to leads a very simple result can be obtained for the 
average current: 

^-^'^'^^ ^ {'^i I ^(I-A(e,t))[/°(6) I . (19) 

In the noninteracting single-level case (IniA(e, t)) is given by 
n A ( A\ r ^ J,^[(Ao-A^/^.)/a;] 

We shall use these expressions below. 



4.4 Resonant tunneling with phonons 

Resonant tunneling diodes (see Fig. |^) are important both for technical ap- 
plications, such as oscillators or infra-red detectors, but they are equally im- 
portant conceptually as a very clear-cut example of a system exhibiting in- 
teracting, far-from equilibrium transport phenomena. In this subsection we 
consider electron-phonon interactions. As explained in the Introduction, the 
IV-characteristic is dominated by a maximum when the quasibound state is 
aligned with thf; energies of the incoming electrons. Experimentally, however, 
it was found lB that the current-voltage characteristic could also display a 
secondary maximum, i.e. a satellite of the main feature. The reason is interac- 
tion with phonons: an electron, which approaches the double-barrier structure 
with a nonresonant energy, can be tuned in the energy of the quasibound state 
in the quantum well by emitting (or absorbing) an optical phonon, and thus 
become resonant with an enhanced transmission probability, and increased 
current. The central-region Hamiltonian (^ is a mathematical formulation for 
this physical picture. 

For simplicity, we consider only energy-independent level- widths Tl and 
Fij, when the current (^ becomes 

•^^hWTf^ y J Jte^^'Ait), (21) 

where A{t) = i[G'^{t) — G°'{t)] is the interacting spectral function. In general, 
an exact evaluation of A{t) is not possible, however, if one neglects the coupling 
to the contacts (or treats the coupling phenomenologically by introducing jan 
exponential decay in A{t)), G^{t) [and hence A{t)] can be calculated exactlyE3. 



10 



A very convenient way to calculate A{t) is to use the linked-cluster theorem 
E3. We write the electron-phonon interaction in the interaction picture (for 
simplicity we again consider just one level), 

= dUj^M^ [at e— ^* + a^e^'*] . (22) 
q 

^in-t^^W '^^^ then be viewed as an effective one-electron Hamiltonian in a 
time-dependent Schrodinger equation, the solution of which gives A(t) (the 
averaging is now over the phonon subsystem): 



A(t) =e-"«*(T{exp[-z f dhH^;"'' 

JQ 



itl)]}) 



0* exp[-i dt, 1^ dt,{T{H^^;^\t,)H^^;^\t,)})] , (23) 



where T is the time-ordering operator, and we used the linked-cluster theorem 
on the second line. The time-ordered product is nothing but the free phonon 
Green function: 

where N = l/[exp(?iti;^/3) — 1]. The time-integrals are easily worked out with 
the result 

A{t) = exp[-it{eo - A) - $(t) - r|t|/2] , (25) 



where 



A = 5:-^, (26) 



q 



$(i) = ^ ^[iVq(l - e^"-**) + (TVq + 1)(1 - e-'""*)] . (27) 

q '^1 

When substituted in the expression for current, one recovers the result of 
Wingreen et aleil, which originally was derived by analyzing the much more 
complex two-particle Green function 

G(r, s, t) = e{s)e{t){d{T - s)d\T)d{t)d^O)). (28) 

The advantage of the method presented here is that one only needs the single- 
particle Green function to use the interacting current formula (H). Other 
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systematic approaches to the single-particle Green function can therefore be di- 
rectly applied to the curmnt (e.g., perturbation theory in the tmmcling Hamil- 
tonian; Anda and FloreO, Hyldgaard et alc5) . The model studied in this sec- 
tion is quite flexible and it has been used to describe many other interaction 
mechanisms in double-barrier systems, such as light assisted magnetotunneling 



4-5 Multiterminal generalization and displacement current 

We now return to the issues of current partition and displacement current in a 
multiprobe sample, mentioned in the Introduction. The dc situation has been 
exhaustively analyzed, starting from the seminal work of Biittiker et alc3. The 
time-dependent case is more complicated: not only has the tunneling current 
to be considered, but also the displacement contributions must be accounted 
for. These issues are essential for two reasons: i) Current conservation and 
ii) Gauge invariance. (Gauge invariance means that a uniform shift of all 
voltages should not affect the final results. jL,A linear, low frequency theory 
was developed by Biittiker and co-workers cZI, based on an generalization of 
the scattering matrix formulation of transport theory. Nonequilibrium Green 
functions allow, at least in principle, the analysis of high frequencies and far- 
from equilibrium situations. Here we will briefly review recent progress within 
this formulation. Stafford discussed the current partition in dn. but under 
nonlixLear conditions Cj, and recent work by Anantram and DattaEJ and Wang 
et alEj have led to an ac generalization. Concerning the tunneling contribution, 
it is a trivial matter to generalize from the two-probe geometry to a multiprobe 
system: one merely replaces the indices L, R by a, where a labels the probes, 
and considers the current Ja{t) (or Jaiuj), whichever is more convenient). The 
dynamic conductance Gap due to the tunneling current is defined as 



In the time-dependent case the tunneling currents do not add up to zero, due 
to charge accumulation/depletion. The total current, however, is conserved: 





(29) 




(30) 



where 



/dE 
—TT[g<{E + u:,E)]V0 



(31) 
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is the accumulated charge in the scattering region; g^{E,E') is the double 
Fourier trantifnrrn nf the Small signal component of the full Green function 
G^, See Refs.p-E^l]. The total current in probe a is J^°* = Jq + J^, where 
^^Ja — J'^ — dQ/dt, and current conservation means J2a'^a'^ ~ ^- 
ditional information is required to partition J'^, because only the sum of the 
various displacement currents is known via Eq.(^). In a model where cou- 
pling constants between the central region and the contacts are independent of 
energy, one can readily do this partitioning: = i^a/ ^/3^i3)J'^- In a more 
elaborate model Wang et alE^I have outlined a procedure how this partitioning 
can be carried out; the analysis is based on requirements of charge conserva- 
tion, G^°^ = 0, and gauge invariance, G^°^ = 0. The end result for the 
dynamical conductance is 



'-'a/Sy^) - ~ ^f3 Qd 



(32) 



where 



Gd 
f3 



-quj 



(33) 



The repilt (32) formally agrees with the scattering matrix results of Biittiker 
et al. e3, but now the various terms are expressed in terms of nonequilib- 
rium Green functions, and hence powerful techniques for evaluating them are 
available. Clearly the final word is not said in this rapidly developing subfield. 



^.6 Phase- measurement of a quantum dot 

Measuring the phase of a transmission coefficient in contrast to the amplitiidfi 
(which determines conductance) has become possible only quite recently Eilo. 
The experimental protocol can be summarized as follows: A magneto-transport 
measurement is performed on an Aharonov-Bohm ring with a quantum dot 
fabricated in one of its arms. If the quantum dot supports coherent transport, 
the transmission amplitudes through the two arms interfere. A magnetic field 
induces a relative phase change, 27r<i>/<i>o, between the two transmission am- 
plitudes, io and tqo, leading to an oscillatory component to the conductance 
g[B) = {e^/h)TiB), with 

TiB) ^ r(°) + 2Re{t5tQ„e2"*/*"} + (34) 

where $ is the flux threading the ring, $o hc/e is the flux quantum, and 
where the dots represent higher harmonics due to multiple reflections. The 
amplitudes to and {qd give the coherent parts of the two sets of paths joining 
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Figure 5: Schematic lay-out of the proposed multi- terminal double-slit interference experi- 
ment. The device consists of an Aharonov-Bohm ring with a quantum dot (QD) in one arm, 
defined with metallic gates (shaded areas) on the 2D electron gas. Electron paths (shown 
as dashed lines) originating from the emitter (E) interfere at the collector (C). A time- 
dependent voltage V{t) is applied to the quantum dot via a side gate. The reflector gates 
R1...4, shown as white areas, direct multiply reflected paths to the common base (B). thus 
preventing them from contributing to the interference signal, in accordance with Eq. (M). 



the emitter and the collector; the incoherent components lead to a structureless 
background signal, which can be neglected in the forthcoming analysis. In the 
experiments, an oscillatory component in magnetoconductance of this form was 
clearly obsericed-thus demonstrating coherent transnaission through the arm 
with the dot E]e3. In the experiment of Yacoby et al. Lil, the Aharonov-Bohm 
phase could take on only two values, and tt, as a consequence of microre- 
pzErsibility in a two-terminal geometry. The second generation of experiments 
E3, in a four-terminal geometry, allowed the determination of the continu- 
ous phase shift of the transmission amplitude through the dot. The success 
of these experiments suggests applications to other phase-coherent transport 
processes. One particular example which has been of considerable recent in- 
terest is photon-assisted tunneling. While photon-assisted tunneling (PAT) is 
intrinsically a coherent phenomenon, existing measurements of PAT are insen- 
sitive to the phase of the transmitted electrons and do not directly demonstrate 
coherence in the presence oLthe time-dependent field. Here we give a brief ac- 
count of a recent proposalEj for a measurement of photon-assisted tunneling 
through a quantum dot in the mesoscopic double-slit geometry described above 
(see Fig. |^). The prpppsed experiment is a combination of the experiments of 
Kouwenhoven et al. eSeS where a microwave modulated side-gate voltage gave 
rise to photon-assisted-tunnieling through a quantum dot, and the interference 
experiments of Refs. andES. 



We can use the results presented in Section 4.2 to construct an expres 



sion for the interference signal. We focus on transport in the neighborhood 
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of a single Coulomb oscillation peak associated with a single nondegenerate 
electronic level of the quantum dot. The effect of the ac side-gate voltage is 
described entirely through the time-dependent energy of this level, which has 
the familiar form 

eit) = eoiVs) + V^c cos Lot, (35) 

where the static energy of the level eo depends on the dc side-gate voltage Vs. 
All other levels on the dot can be neglected provided the ac amplitude, Vnc, 
and the photon energy, huj, are small compared to the level spacing on the dot. 

The energy-dependence of the coherent part of the transmission ampli- 
tude Iqo (e) through the arm containing the quantum dot is determined by the 
transmission amplitude tQD(e) through the dot, tQD(e) oc tQo[e). In the absence 
of an ac potential, a suitable model for the dot transmission amplitude is the 
Breit-Wigner form, 

where F = Fl -I- Fi? is the full width at half maximum of the resonance on the 
dot due to tunneling to the left and right leads. Eq. (|6|) implies a contin- 
uous phase accumulation of tt in the transmission amplitude as the Coulomb 
blockade peak is traversed. (Note that the Breit-Wigner form is exact for a 
noninteracting system with F independent of energy.) 

In the dynamic case, the simple Breit-Wigner description- must be gen- 
eralized, and the object to evaluate is the ^-Matrix element BcJ. Provided 
interactions in the leads can be neglected, the elastic transmission amplitude 
tooie) can be written as the energy conserving part of the S'-Matrix between 
the left lead and the right lead 

lim (e', i?|5|e, L) = <5(e' - e)tQ^{t), (37) 

where 

tQn{^)=-i^TLVR{A{e,t)), (38) 
with (compare to Sect. [4.3| ) 

A;— — OO 

The result for the transmission amphtudc at finite temperatures is 



k— — oo 



AttT 



2 2'KkBT 



(40) 
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Figure 6: Temperature dependence of the phase shift A0 (top panel) and the square of 
the amphtude (bottom) of Iqo- The level-width is r/2 = 0.1, in terms of which the other 
parameters are Vkc = 1-0, u) = 1.0, and T = (solid line), 0.1 (dashed line), 0.5 (dash-dotted 
line). For comparison, the T = time-independent results are shown as dots. 



where tp' is the derivative of the digamma function, and /i is the chemical 
potential in the leads. 

Figure shows the computed magnitude of Iqc (bottom) and its phase 
(top), as a function of the level energy eo(Ki). As compared to the time- 
independent case (shown as a dotted line), several features are noteworthy. 
The magnitude of Iqo shows photonic side-bands, reminisceat of those seen in 
transmission through a microwave modulated quantum dotEj. However, there 
is an important difference from the usual case of photon-assisted tunneling. 
The amplitude of the Aharonov-Bohm oscillation is sensitive only to the time 
average of the transmission amplitude Iqo- Hence only elastic transmission 
through the dot contributes, i.e., the net number of photons absorbed from 
the ac field must be zero. The sideband at say e — eo{Vs) — fiu corresponds 
to a process in which an electron first absorbs a photon to become resonant 
at energy eo{Vs)^ and subsequently reemits the photon to return to its original 
energy. Perhaps most interesting are the features appearing in the phase: the 
phase shift shows a non-monotonic behavior, with pronounced resonances lo- 
cated at the energies corresponding to the photonic side-bands. The strengths 
of these phase resonances are strongly dependent on the ac amplitude T4c, 
and in Fig. ^ we highlight an interesting consequence of Eq. (^0|) : it is possible 
to entirely quench the main transmission peak (bottom panel), or change the 
sign of the slope of the phase at resonance by adjusting the ratio Vac/huj to 
coincide with a zero of the Bessel function Jq (top) . This phenomenon is math- 
ematically analogous to the recjeiitly observed absolute negative conductivity 
in THz-irradiated superlatticesE3; in our case, however, it is the phase rather 
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Figure 7: Temperature dependence of the phase shift (top panel) and the amphtude (bottom) 
oitQD, for Vac = 2.405, ui = 1.0, V/2 = 0.1, T = 0,0.1,0.5. 



than the current that displays this behavior. For a reader interested in further 




dopments along these lines, we direct attention to two very recent papers 



4.7 Time- dependent Kondo physics 

The Kondo effect is undoubtedly one of the most studied and best understood 
problems of many-body physics. Initially, the theory was developed to ex- 
plain the increase_Df resistivity of a bulk metal with magnetic impurities at 
low temperatures E3. More recpnJJy, the observation of Kondo effect in several 
experiments on quantum dot£3l£j has demonstrated that these systems can 
serve as an important new tool in the study of strongly correlated systems. 
Unlike magnetic impurities in metals, the physical parameters of the quantum 
dot can be varied continuously, which allows, for example, systematic experi- 
mental study of the crossover between the Kondo, the mixed- valence, and the 
non-Kondo regimes. These aforementioned. C|Xp|priments were in good agree- 
ment with earlier theoretical predictions e3'E3E3|£3. It is therefore natural to 
ask: can Kondo physics,pbe probed in time-dependent experiments? Already 
before th r.rxnqprnentsE^i^j a few theoretical papers addressed ac-driven Kondo 
systemsE300E3, but thij-Last-two years have witnessed a flurry of theoretical 
activity, see, e.g., Refs.pOBB]. 

A full discussion of these papers would merit a review of its own. Our 
present goal is however much more modest. This is a rapidly developing area 
and still ridden with some controversies, which will provide fuel for continued 
intense research efforts. Important for this volume is that the main technical 
work-horse has been the nonequilibrium Green function technique. It is there- 
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fore hoped that the reader will get a feeling of the present excitement from 
the following brief discussion on some of the recent achievements and open 
problems related to the ac-Kondo problem. 

The static Anderson model is solvable with Bethe AnsatzlIS or by quantum 
Monte Carlo methodsllj, but a reliable and simple method to obtain dynamical 
properties at low temperatures in the whole range of interactions (U/T) is 
not available. A number of different methods have been used to attack the 
problem. For example, the non-crossing method (NCA) may be used following 
an exact trajasformation of the U = oo Anderson model into a slave-boson 
HamiltonianL3. The latter is then solved self-consistently to second order in the 
tunneling matrix elements Vk ■ The NCA gives reliable results for temperatures 
down to r < Tx, and its time-dependent formulation has been very usefuLin 
applications to problems in surface physicsOLll. This is the approach of RefQ, 
which studies the temporal response of a quantum dot which is suddenly shifted 
into the Kondo regime by a change of a voltage on a nearby gate. Thus, the 
NCA is applied to a calculation of the nonequilibrium spectral density. It 
is suggested that subjecting the quantum dot to a sequence of pulses, and 
analyzing the resulting current as a function of the duration of a pulse, will 
open a window to the build-up mechanism of the many-body correlations, 
responsible for the Kondo effect. 

NCA, however, does not give exact results as T — > 0, and other methods 
are called for. Finite U perturbation theory works for the symmetric case, 
but is known to exhibit anomalies away from this special case. An attempt 
to circumvent these problems is to construct an effective self-energy, which 
smoothly interpolates between known limits U/T and T/U 0, and 

Sua may eliminate some of the problems related to perturbative approaches 
Ea. W|hile a plausible process, interpolation is not rigorous, and the authors 
of RcfEa, suggest a nonequilibrium generalization of the Friedel sum rule to 
construct a consistency check in the computation of the nonequilibrium Green 
functions. The resulting time- averaged density-of-states exhibits a rich struc- 
ture: one may find replicas of the Kondo peak and/or mean-field peaks. Also, 
the conductivity is found to be strongly affected by the ac-drivinErnand, in- 
terestingly, that it cannot be described by the Tien-Gordon modelEJ, i.e., by 
single-particle ac-assisted tunneling. The observability of the satellites-pf the 
Kondo peak is, however, brought into question by Kaminski et al. Ell, who 
point out the that the coherence of the many-body correlations responsible 
for the Kondo anomaly is fragile against a spin-Mip cotunneling process, which 
may occur already at relatively low frequencies Ed. It is clear that experiments 
exploring the ac-Kondo system would be extremely welcome! e3 
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4-8 Semiconductor- superconductor hybrid systems 

Materials technology has in recent years advanced to a stage where high qual- 
ity hybrid structures can be fabricated. By hybrid structures we understand 
combinations of normal metals (N) and superconductors (S) together with 
mesoscopic structures, such as quantum dots (QD). Thus, in Fig. |^ cither of 
the contacts, or the quantum dot itself could be superconducting. Obviously 
there are many possible combinations and here we just dismss one of them, 
namely a N-QD-S structure, recently studied by Sun et alEda. 

ae superconducting contact allows the possibility of an Andreev reflec- 
an electron approaching the superconductor with subgap energy may 
be transferred into the condensate with a simultaneous creation of a back- 
propagating hole. This process reflectsJiself in a number of ways, for exam- 
ple as bound states in a S-N-S system e£|, or the even-odd parity asymmetry 
and the CaLLlomhJblockade of the Andreev reflection in S-SQD-S or N-SQD- 
N systems eZIO'E3l3. Traditionally rf radiation has been one of the standard 
probes for studying superconductivity, and a large number of experimental 
results have been reported on "simple" hybird systems. Here we discuss the 
case where a single level (with spin) couples to both radiation and N and S 
contacts a. (It should be noted that a strongly inteeacting QD in this con- 
figuration would bring us again to the Kondo realm Ell.) The Hamiltonian is 
thus H = HL{t) + i?ccn(t) + Hr + Ht,l + Ht,r, where Hl, H^en and Ht,l 
are standard but the presence of the superconductor requires the following 
modification for the terms referring to the right-hand side of the system: 

Hr = J2 ^v'^lo^Po + [^*CpiC-P>T + ^cLp^^c], J (41) 

per p 

ffr.fl = E [^fee-^«'ct^d. + y;e— ^«*4cp.] , (42) 



pa 



where the voltage of th£_^uperconducting contact Vr appears in the phase of 
the tunneling coupling l£3. This can be viewed as a mean field Hamiltonian 
leading to the Bogoliubov-de Gennes equations. 

The derivation of the time-dependent current proceeds along the same 
lines as in Sect. ||. In particular, the central result Eq.([Tl|) is still valid, if one 
interprets the central region Green functions as Nambu matrices, 

^ (*'^)-^l,(cj(i')cl(t))(c,(Octi(t)j ' 

and analogous definitions hold for the retarded Green functions and self-energy 
functionals. The ensuing calculations are quite complicated, even in the case 
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of a noninteracting quantum dotEj, and we focus on some of the main physical 
conclusions. First of all, one must distinguish on what part of the system 
the time-variation is affecting, and whether the photon energy is bigger or 
smaller than the energy gap. As an example, let us consider the case when 
huj < A, and that the external radiation affects the superconductor (by gauge 
invariance, this is equivalent to the case when the normal contact and the dot 
are affected, but the superconductor is not). The dominating contribution 
to the current arises from a photon-assisted Andreev term (PA AT), which 
has severiai interesting properties. Among these is a possibility of electron 
pumping l3: the electrons move from a lower potential to higher potential 
by absorbing-4)hotons, leading to a negative current (or, absolute negative 
conductivity eSI) . If, on the other hand, h > A, the normal PAT-processes 
contribute substantially. In general, much interesting structure is seen as a 
function of the gate voltage, which can be used to tune the energy levels in 
the quantum dot—and the structure is quite different from the one seen in 
N-QD-N systemM 

5 Concluding remarks 

In this review we have attempted to give some insight to recent developments 
in time-dependent transport in mesoscopic systems, treated by nonequilibrium 
Green function methods. The review is by no means exhaustive: maav-im- 
portant and interesting problems have pot been covered, such as noise or 
surface acoustic wave driven transport e3 (and there are many other topics as 
well) . The scope of the review has not allowed an in-depth treatment of many 
of the topics, but the author's hope is that the incomplete treatment has raised 
the reader's curiosity, and perhaps thereby attracted new researchers into this 
dynamic field of research. 
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